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g\ ■ Abstract 

a^ 

K^i We study color confinement properties of the niulti-instanton system, wliich 

C^ I seems to carry an essence of the nonperturbative QCD vacuum. Here we 

assume that the multi-instanton system is characterized by the infrared sup- 
pression of instantons as f{p) -^ p^^ for large size p. We first investigate 
a monopole-clustering appearing in the maximally abelian (MA) gauge by 
'^ ! considering the correspondence between instantons and monopoles. In or- 

ly^ I der to clarify the infrared monopole properties, we make the "block-spin" 

^D ' transformation for monopole currents. The feature of monopole trajectories 

changes drastically with the instanton density. At a high instanton density. 



^ 



o 

^^ . there appears one very long and highly complicated monopole loop covering 

2^ I the entire physical vacuum. Such a global network of long-monopole loops 

resembles the lattice QCD result in the MA gauge. Second, we observe that 
^ ■ the SU(2) Wilson loop obeys an area law and the static quark potential is 

-L . approximately proportional to the distance R between quark and anti-quark 

D I in the multi-instanton system using the SU(2) lattice with a total volume of 

'^ ' V = (lOfm)^ and a lattice spacing of a = 0.05fm. We extract the string 

tension from the 5 x 10^ measurements of Wilson loops. With an instanton 

^ ■ density of (N/V) = (1/fm)^ and a average instanton size of p = 0.4 fm, the 

H , multi-instanton system provides the string tension of about 0.4 GeV/fm. 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) has been established as the fundamental theory of the 
strong interaction. In the infrared region, there appear various nonperturbative phenomena 
such as color confinement and dynamical chiral-symmetry breaking. Since the QCD vacuum 
is composed of gluon fields interacting in a highly complicated way, it is hard to understand 
these phenomena from perturbative points of view. On the other hand, a topological aspect 
may provide a useful approach for descriptions of the QCD vacuum. Actually, there appear 
two non-trivial topological objects, instantons and monopoles, due to the non-linearity of 
QCD. 

The instanton configuration discovered in 1975 by Belavin, Polyakov, Shvarts and Tyup- 
kin in is a classical and self-dual solution of the Euclidean field equation in the Yang-Mills 
theory. The appearance of instantons corresponds to a homotopy group, n3(SU(A^c)) = Zoo 
0. Instantons are important for nonperturbative phenomena related to the Ua(1) anomaly 
and the large rj' mass P]. Chiral-symmetry breaking could also be interpreted as the in- 
stanton effect [^|-|§]. However, until now there has been no evidence that instantons have 
anything to do with color confinement in the 4-dimensional gauge theory, although Polyakov 
discovered that instantons cause confinement in certain 3-dimensional Georgi-Glashow mod- 
els 0. With recent computational progress, instanton properties are investigated by lattice 
QCD simulations based on the cooling procedures, which are achieved by an artificial reduc- 
tion of the local lattice action. This method allows one to eliminate short-range quantum 
fluctuations of gluon fields and to extract only topological excitations like instantons from 
the nonperturbative QCD vacuum [p|-p!3|]. These investigations provide us the average in- 



stanton size p ^ (0.33 — 0.4) fni and the instanton number density (N/V) ~ 1 f m '^ IH"14 



Although the existence of instantons with high density is agreed, the detailed numbers for 
p and the size distribution depend largely on the cooling procedure. These scale parameters 
suggest that the QCD vacuum is the dense matter of instantons and anti-instantons. 

In 1981, 't Hooft proposed the abelian gauge, where the color-magnetic monopole appears 
as a relevant degree of freedom for the description of color confinement |15| . In the abelian 



gauge, the SU(iVc) nonabelian gauge theory is reduced to the 11(1)^'="^ abelian gauge theory 
with color-magnetic monopoles. The appearance of magnetic monopoles corresponds to 
another homotopy group, n2(SU(A^c)/U(l)^''~^) = 2^^="^, which is different from that of 
instantons. In the abelian gauge the color confinement mechanism can be interpreted as 
the dual Meissner effect due to monopole condensation, which is a dual version of Cooper 
pair condensation in the ordinary superconductivity. Such a dual superconductor picture 
for color confinement was proposed by Nambu, 't Hooft and Mandelstam in the middle of 



1970's ||16 HTq] - By lattice QCD simulations, it is observed that large monopole clustering 
covers the entire physical vacuum in the confinement phase, which is identified as a signal 
of monopole condensation being responsible for confinement |T^J2^. Many studies indicate 



that the monopole is a relevant degree of freedom for color confinement and chiral-symmetry 
breaking |l]-g|]. 

Recent studies show remarkable facts that instantons are directly related to monopoles 
|p9|^7| in the abelian gauge, although these topological objects belong to different homo- 
topy group. For example, in the Polyakov-like gauge [^], monopoles appear due to the 
existence of the hedgehog configurations near instanton centers. Such correlations indicate 
that instantons may be important for the promotion of long monopole loops. Therefore, we 



take the inulti-instanton system instead of the QCD vacuum in order to clarify monopole 
clustering and confinement properties in terms of instantons. 

II. TOPOLOGICAL OBJECTS IN THE QCD VACUUM 

The instanton is a classical and non-trivial solution of the Euclidean Yang-Mills theory, 
whose action is written as 5* = / d'^xG'^^GI^^^/Ag'^. Here, in order to make the notation and 
the discussion simpler, we take the SU(2) case and represent A^{x) = igA'^{x)T°-/2 and 
G^y{x) = ■i^'G" (x)r'^/2, which are defined as the anti-hermite variables. Instanton and 
anti- instanton configurations are characterized by the (anti-) self-duality condition, 

by using the dual field strength G^,y = ^s^^a/sGais- This condition provides the minimal 
action S = ^ \ Q \. Here, the topological charge is defined as Q = /d^x{Gj^^G^^}/327r^ 
These instanton solutions satisfy the general Yang-Mills field equations automatically, 
D^G^y = ^D^G^y = 0. The self-dual solution with Q = 1 in the singular gauge [Q is 
written as 

Here, the instanton solutions have several collective modes related to the size p and the 
position z of instanton, which are described by one and four parameters, respectively. For 
pure SU(2) gauge theory, the instanton solution can be rotated in color space by the color 
orientation matrix O, which is characterized by 3 parameters as the Euler angle. The 't Hooft 
symbol ff^^ is defined as 

-b^!/ _ _-feiy^ _ I ^ 1 /i, z/ = l,2,3 ,s 

^ - n -\-5^^ z/ = 4. ^^^ 

The anti-self-dual solution A^ is obtained by replacing ff^^ to rf'^^ = {"^y '^^'^ V^^'^ i^i 
Eq.(i. 

Now, let us discuss the appearance of QCD-monopoles in the abelian gauge. We consider 
the maximally abelian (MA) gauge |T^, which is defined in the Euclidean SU(2) QCD by 
minimizing the variable 

R.n[A,]= I d'x{{Al{x)f + {Al{x)f} (4) 



under the gauge transformation VL{x). In the MA gauge, the off-diagonal fields are suppressed 
by the gauge transformation and the full gauge field A^ = A'^^^ behaves as the abelian 
gauge field a^ = A^^^ approximately. In this paper, we adopt the MA gauge, where the 
abelian dominance holds for the nonperturbative QCD phenomena [^]. Under the gauge 
transformation, the gauge field transforms as 

A^{x) ^ A]^{x) = n{x)iA^ix) + d^MxY. (5) 



After the gauge fixing, there only remains the abelian gauge symmetry 11(1)'^'="^ C SU(A''c). 
When we meet hedgehog configurations like instantons, the multi-valuedness happens to the 
gauge function Q{x) and the gauge transformation of Eq.(|^) develops a singularity. This 
singularity leads to the monopole current fc^. The field strength is defined generally as 
Gfj,i, = [Dfj,,D,,] — [d^,d^] [Q, which returns to the standard definition G^,^ = [Dfj_,D^] = 
d^Ay — dyAfj, + [Afj,, Ay] for the regular case. Under the singular gauge transformation Q{x), 
the field strength transforms as 

G;^, = nG,^Q^ = d,A^ - d^A^^ + [A% A^] - Q[d„ d^p. (6) 

Since the last term is diagonal, abelian field strength is naturally obtained as 

F^,u = d^ay - dya^ -VL[d^, dy\VL\ (7) 

by performing the abelian projection, A^ — > a^ = tr(74^r^)^. The abelian Bianchi identity 
is broken due to the existence of the last term in Eq. (0) , and the magnetic current fc^ (x) is 
obtained as 

F(x) = -e^uapdyF^p = --e^ya(idyn[da, df3]^K (8) 

The obvious consequence of the monopole current conservation, d^k'^ = 0, means that 
monopole currents form closed loops. 

III. CORRELATION BETWEEN INSTANTON AND MONOPOLE FOR COLOR 

CONFINEMENT 

Recently, both analytical and lattice studies showed a strong correlation between instan- 
tons and monopoles in the abelian projected theory of QCD [p9|-^. Let us briefly review 



recent studies on the correlation between these topological objects in the MA gauge ||35| , |39 
The minimizing condition of i?ch[^^] in Eq.(|) satisfies the local condition, 

(^mT 4)4 = 0, A^^Al±^Al. (9) 

It is noted that self-dual solutions like a instanton satisfy the stationary condition {^ auto- 
matically. However, this stationary condition is not sufficient to realize the MA gauge. In 
the MA gauge the functional i?cfc[^^] must be minimized in addition. The arbitrary gauge 
choice of a single instanton gauge field leads to a different value of Rch[A^], although the 
instanton gauge field satisfies the condition (P) in any gauge [^. For instance, in the singu- 



lar gauge, which has a point singularity at the instanton center, the instanton configuration 
gives a finite value Rch[At\ = 47r^p^ with p being the instanton size. On the other hand, 

in the nonsingular gauge, there appears the divergence as i?c/i[^J]] = 2 J d'^x ,^1 2\2 — *■ oo. 
Therefore, it is necessary to consider the minimization of the functional Rch[A^] for the MA 
gauge fixing. 

Brower et al. consider the single instanton configuration in an "intermediate" gauge 
where the gauge field has the singularity at the closed loop with a radius R around its 
center [B^]. This singularity leads a closed monopole loop of radius R around the instanton. 



The parameter R should be decided by minimizing the functional Rch [^/x] • It is numerically 



shown that the monopole loop with radius R prefers to shrink to a point |^. This fact is 
natural from the following consideration. There is no definite direction both in the single 
instanton configuration and in the MA gauge condition due to the 4-dimensional rotation 
invariance. Then, the normal vector of the monopole loop cannot be fixed and the loop 
ought to shrink to a point. However, the R dependence of Rch[-A.^ is found to be extremely 
small. Therefore, some small perturbations as the presence of the finite volume cause a 
closed monopole loop with a non-zero radius R due to the fragile nature of the monopole 
point solution, which is first reported by Hart and Teper using the lattice simulation []35[ . 



In the QCD vacuum, there are actually many instantons and anti-instantons, and here the 
appearance of monopole loops should be influenced by other instantons and anti-instantons. 
Therefore, we would like to investigate the monopole loops by using a more realistic multi- 
instanton system defined on the R^ space. 

We first study the monopole clustering as a signal of monopole condensation in terms 
of instantons, which is motivated by these strong correlations between instantons and 
monopoles in the MA gauge. Second, we calculate the static quark potential in the multi- 
instanton system by using the Wilson loop in order to clarify the relation of instantons with 
color confinement. 



IV. MULTI-INSTANTON MODEL 

In this section, we would like to model the nonperturbative QCD vacuum in terms of 
instantons. The QCD vacuum possesses the gluon condensate [^, which relates the number 
of instantons and anti-instantons |^. If the total action can be estimated as the sum of 
individual instanton actions, the gluon condensate is proportional to the instanton density 
as (G'^^G'JJ^)/327r^ ~ {NjV). The QCD sum rule provides the phenomenological value 
(G'^j,G^j,)/327r^ ^ (200 MeV)"^, and hence the average instanton density becomes [N/V] ~ 
(1/fm)'^ by the above assumption. The lattice QCD simulations suggest that the QCD 
vacuum is saturated with many instantons and anti-instantons. Hence, we model the QCD 
vacuum by the multi-instanton ensemble. First of all, we consider the partition function 
^Inst of the single instanton as the basic ingredient of the multi-instanton theory. Using the 
collective coordinates, p, z^, O, the single-instanton partition function is expressed as 



C{N,)^ Stt^ ,2^^,,^^,, _, Stt^ 



with b = yA"c- Here, /o(p) is the single-instanton weight function in the one-loop approx- 
imation [^,^]. This weight function has a scale M corresponding to the scale invariance 
breaking by the trace anomaly in QCD. The bare coupling g'^{M) is given at this scale M. 
The weight function /o(p) with Nc> 2 increases with the instanton size p, and the infrared 
divergence appears in Zl^^-f. within the one-loop approximation. 

We consider now the multi-instanton system. Based on the single-instanton partition 
function, the multi-instanton partition function is expressed as 



N++N. 

n 

N+N. ^» + . ^»-. „^| 



1 1 ^++^- f f f 

^ ^ Y^ AMAM n J d^Zn J dpnj dOnfoiPn)exp{-UintiZn,pn,On)}, (H) 



where Uint denotes the interaction between instantons (anti-instantons) and depends gener- 
ally on Zn, Pn, On- Here, we follow the standard method, where the interaction Uint is taken 



to depend only on the size p ||5^-|54|. Therefore, we consider now f{p) = fo{p)exp{—Uint} 
as the instanton size distribution. This interaction is found to be repulsive, which sup- 
presses the appearance of large size instantons |B2|-E^. We regard the positions and the 



color orientations as random variables in our calculation. 

For the small instanton, the perturbative scheme is valid and the interaction Uint can be 
neglected. Hence, the instanton size distribution f{p) behaves as the single instanton weight 
/o(p), 

/(p)^ const ■/-^ (12) 

On the other hand, for the large instanton, the direct estimation of f{p) is very complicated 
due to the nonperturbative properties. The analytical studies and the numerical lattice 
QCD calculations fS^JS^ suggest a strong suppression of the large size instanton as 

fip)'^ const -p-^ (13) 

which is caused by the repulsive force in the infrared region. The ordinary instanton liquid 
model suggests z/ = 5 p6| . 

To connect the two tendencies in Eqs.(0) and (|13|) smoothly, we take the size distribution 

as 

where pi denotes the infrared size parameter and p2 the ultraviolet size parameter. These 



two parameters are fixed by the average instanton size p = J^ dppf(p) = 0.4 fm |^^ and 
the normalization condition /q°° dpf{p) = 1. 

The above discussion leads to the use of the partition, function 

1 1 ^++^- f f°° f 1 



As for the gluon field A^ of the multi-instanton system, we take the sum ansatz 



\i^) = Yl ^i(^5 Zk, pk, Ok) + Yl ^i(^5 z-k, pk, Ok), (16) 

k k 

which is constructed by the instanton and anti-instanton solutions in the singular gauge. 
Using the variational treatment with this ansatz, Diakonov and Petrov have shown the 
appearance of the repulsive force between instanton and anti-instanton in the infrared region 
||. 

In actual calculations, we generate the ensemble of instantons and anti-instantons with 
random centers Zk on the 4-dimensional Euclidean continuum space. The color orienta- 
tions Ok are taken randomly. The instanton sizes pk are randomly chosen following the 
size distribution /(p) in Eq.(|l^). In the simple sum ansatz, the interaction among these 
pseudoparticles is supposed to be included effectively in the instanton size distribution /(p). 
The instanton size distribution for the cases of z/ = 5 and 3 are shown in Fig.|l]. 

6 



V. MONOPOLE CLUSTERING IN MULTI-INSTANTON SYSTEM 

Based on the lattice gauge theory, we investigate monopole-loop distributions induced by 
instantons after the MA gauge fixing [|^,Q. We introduce a lattice on the multi-instanton 



configuration and define the link variable ?7^(s) = exp(zayl^(s)), where A^{s) is provided on 
each link from Eg. (]TB|) . In the actual calculation, we use the 32^ lattice with the lattice- 
spacing of a = 0.125 fm. 

First of all, let us discuss the boundary condition of the multi-instanton configuration. 
If instantons and anti-instantons are generated only in the finite volume V = (32 x a)*^ = 
(4.0 fm)^, the strength of the gluon field is insufficient near the border of the volume V. 
To avoid this border problem, we adopt the periodic boundary condition, which ensures 
adequate contributions also near the border. Here, instantons are assumed to appear with 
the periodic position, size and color orientation out of the finite volume V. A schematic view 
of this periodic condition of an instanton ensemble sliced into 2-dimensional plane is shown 
in Fig.^]. This means the consideration of the instantons in the 3^ boxes in the 4-dimensional 
case: it is 3^ = 9 boxes in the 2-dimensional case as shown in Fig.|^. In actual calculation, 
we construct link variables by considering all the instantons in the neighboring boxes. 

Now, we apply the MA gauge fixing [|19[ by maximizing 



R = T.^r[U,{syul{sy] = 2E[1 - HKisyf + (Ulis))'}], (17) 



S,fJ. 



with f/^(s) = Ujl{s) +iT^'W{s). The maximization of R corresponds to the lattice expression 
of the minimization condition of Rch[Afj] in Eq.(^). In the MA gauge, the SU(2) link variable 
U^{s) is decomposed as 

U,{s) = M,{s)u,{s) = ( V ' ' ;-- ' ; ^ 77,,, I r Q ^_i(.) ) , (18) 




where the abelian angle variable 6'^(s) and the non-abelian variable c^(s) are defined in 



terms of f/^(s) as 

tane,(s) = -^, c,(.)e^^^W = [-U^is) + zUl{s)]. (19) 

It is obvious from the expression of Eq.(^) that the off-diagonal parts t/^(s) and f/^(s) 
of gluon fields are minimized by the MA gauge transformation. Therefore, full SU(2) link 
variables would be approximated as U(l) link variables, f/^(s) ~ ti^(s), in the MA gauge. 



Monopole currents can be defined by using u^(s) following DeGrand and Toussaint [|59 
Using a forward derivative 9^/(s) = f{s + ft) — f{s) with unit vector ft, the 2-form of the 
lattice formulation, 0^y{s) = d^9y{s) — duO^{s), is decomposed as 

9^,{s)=9^,{s) + 27rn^,{s) (20) 

with 0^y{s) = mod2n0^iu ^ ("^r, tt] and n^y{s) G Z. Here, 9^u{s) and 2-nn^i,{s) correspond to 
the regular field strength and the singular Dirac string part, respectively. Since the abelian 
Bianchi identity is broken, the monopole current k^^{*s) can be defined on the dual link (*s, /i) 
as 



k^,{*s) = —e^,uai3du0af3{.s + jj) = -dyh^^{*s) (21) 

where h^u{s) = \£ ^luapfT'api.s + fi). The obvious current-conservation law d'^k^{*s) = leads 
to the closed monopole loop in the 4- dimensional space. Here, d' denotes a backward 
derivative. 

Since we are interested in the infrared behavior of monopoles in multi-instanton config- 
urations, it is useful to execute a "block-spin" transformation on the dual lattice with the 
scale factor A/". The "block-spin" transformation removes small monopole loops as numerical 
noises and keeps its global structure. Here, the A/"^ extended monopole is defined as 

n— 1 

naisis) = J2 Raisins + ip + ja) (22) 

i,j=0 



on a sublattice with the spacing of 6 = Afa pDI. Then, the extended monopole current is 
defined as kf/'{*s) = -\e ^^apdun^pis + A)- 

We start with a very dilute instanton system in order to clarify a local correlation between 
instantons and monopoles in the MA gauge. We take a density {N/V) <^ (0.5/fm)^ = 
(100 MeV)^, where instantons are separated completely from each other. From the analytical 
consideration, one may expect that a monopole loop prefers to shrink to a instanton center. 
However, we observe on the lattice that each finite-size monopole loop is localized around 
each instanton center as shown in Fig.^(a). Such a small monopole but finite-size loop is 
caused by finite lattice spacing and the boundary effect, which is numerically washed out 
by several numbers of "block-spin" transformation. 

In the actual QCD vacuum, instantons saturate the 4-dimensional space and stay close 
to each other with a instanton density of {N/V) = (1/fm)^ [^. Therefore, we consider three 
instanton number-density cases, (N/V)^ = 0.75, 1.00 and 1.25 fm~^, and investigate the 
behaviors of monopole loops. Figure ^a) shows the histograms of unblocked monopole loops. 
The histograms of blocked configurations are shown in Figs, ^(b)-(d), which corresponds to 
the extended monopoles on the sublattice with b = 2a, b = 4a and b = 8a, respectively. 

At the low instanton density where each instanton is isolated, the monopole loop induced 
by the instanton prefers to be localized around each instanton center as shown in Fig.^(a). 
Thus, there appear only relatively short monopole loops. This situation provides a peak at 
zero monopole-loop length in the histogram as shown in Fig.|Kd-i). As the instanton density 
increases, some monopole trajectories tend to hop from one instanton to another nearby 
instanton as shown in Fig.|^(b), and there appear long monopole-loops. Here, a clustering of 
long monopole loops appears and grows gradually to be separated from the small monopole- 
loop part in the histogram. Furthermore, the "block-spin" transformation is available to 
clarify monopole loop behaviors at the large scale by comparing Figs.^(a-iii), (b-iii), (c-iii) 
and (d-iii). Because this procedure removes small-size monopole loops and combine several 
long monopole-loops into one longer loop. At a high instanton density, there appears one 
very long and highly complicated monopole loop in each gauge configuration as shown in 
Fig.^(d)-(iii). Such an appearance of the monopole clustering over the entire physical volume 
can be interpreted as the Kosterlitz-Thouless-type phase transition [p] ,|6^. 



In addition, we discuss the case of the instanton size distribution z/ = 3 with the average 
size p = 0.4 fm as shown in Fig.^. Diakonov pointed out that if the instanton size distribution 
falls off as /(p) ~ 1/p^ in the infrared region, one gets a linear confinement potential []55| , 
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which is based on the formula for the inter-quark potential in Ref. [Q. Since instantons 



tend to have large overlapping from each other for the z/ = 3 case, the simple treatment 
with the sum ansatz, which is based on the assumption of the statistical independence, may 
not be applicable for the construction of the multi-instanton system. However, the infrared 
behavior of the instanton size distribution is interesting for the confinement properties in 
the long-range region. Therefore, we demonstrate also the u = 3 case. The histograms of 
the 1/ = 3 case are qualitatively similar to the corresponding histograms of the z/ = 5 case, 
although there is a quantitative difference between these two cases. 

Lattice QCD simulations at finite temperature suggest that the instanton density is 
largely reduced as the temperature increases. Now, we compare the histograms of the 
monopole loop length in the multi-instanton system with those of the SU(2) lattice QCD 
with 16^ X 4 at different temperatures {f3 = 2.2 and 2.35) in Fig.^ |^1[. The monopole- 
loop distribution at high instanton-density case as shown in Fig.|^(d-iii) resembles with the 
result of the confinement phase (/3 = 2.2) in Fig.|^(b), where many instantons and anti- 
instantons saturate with the instanton density {N/V) ~ (1/fm)^. On the other hand, the 
dilute instanton system as shown in Fig. ^(d-i) is similar to the result in the deconfinement 
phase in Fig. ^(b) obtained by the lattice QCD {(3 = 2.35), where long monopole loops 
disappear. This resemblance seems to indicate that the instanton plays a relevant role on 
the promotion of monopole loops in the confinement phase. 

VI. COLOR CONFINEMENT IN MULTI-INSTANTON SYSTEM 

In the previous section, we have discussed that a high-density instanton system provides 
highly complicated monopole loops. We shall further work out the confining property in the 
multi-instanton system without referring to monopole configurations in the MA gauge [^ . 



For the investigation of confinement properties, it is reasonable to consider the Wilson loop 
with a contour C, which is defined as 

W[C] = trPexpf? <f A^dx^] = trj] U,,is). (23) 

Jc (J 

If the contour C is a rectangle of the dimension T by i?, the Wilson loop is related to the 
static quark potential [^ as 

V{R) = - lim {In W{R, T)/T}. (24) 

To extract the potential without the contamination of excited states, the time T has to be 
large enough compared to the distance R between quark and anti-quark, T ^ R |Q. Hence, 



it is desired to take a large enough T in the calculation of Vr[C]. In the periodic boundary 
conditions, the time distance is limited as a half of Nt x a, where Nt is the number of lattice 
in the time direction and a is the lattice unit. Therefore, we introduce the 36^ lattice in 
the center of a 4-dimensional instanton configurations without the periodic condition for 
the sizes, positions and color orientations of instantons. Here, the large Wilson loop can be 
calculated as T ~ 36 x a because of the non-periodic boundary. A schematic view of this 
condition sliced into 2-dimensional plane is shown in Fig.^. 

In our actual calculation, we consider the instanton configuration with the entire volume 
of V^ = (10 fm)^. We have to fix now the lattice-spacing a. For instance, the instanton density 

9 



{N/V) = (1/fm)^ corresponds the instanton number A^ = A^+ + A^_ = 10^ in the volume 
(10 fni)"^. Since the gluon gauge field is constructed from the summation of instantons and 
anti-instantons as Afj_{s) = J2n+,n^{Ai + ^l); a fine lattice-spacing is necessary to justify 
the continuity of the link variables U{s) = exp[mA^(s)] in above situation. Therefore, we 
actually take a fine lattice spacing of a = 0.05 fm considering the continuity, aA^{s) <^ 1. 

Recent lattice QCD simulations provide us the average instanton size p and the instanton 
density (N/V). The instanton properties are extracted from the smoothed QCD vacuum 
by using several cooling methods. DeGrand et al. investigate a smoothing procedure based 



on the renormalization group equation [Q. In this SU(2) lattice calculation, they find 
the average size of instantons about p = 0.2 fm, which is too small a value than that of 
usual instanton liquid model, at a density of about (N/V) = 2 fm~ . On the other hand, 
de Forcrand et al. use a cooling algorithm based on the improved action with scale invariant 
instanton solution and provide that the instanton size distribution is peaked around 0.43 fm 



I4 . These two cooling methods provide different values about the average instanton size. 
As for the instanton size distribution, there are much effort to extract a size distribution 
from the lattice configurations by various smoothing methods |TT|-|T^ . However, there is no 
consensus on the size distribution either. In our calculation, we take the instanton density 
as (N/V) = (1/fm)^ and the average instanton size p = 0.4 fm. These parameters are 



used in the discussion of chiral symmetry breaking [^ and give us the monopole length 



distribution in the range of those obtained by QCD simulation as shown in the previous 



section [0]. For simplicity, we consider the typical case where the instanton number is 
equal to the anti-instanton one, Nj = Nj = N/2. 

Using the large 36^ lattice without periodic boundary condition, we can make as many 
as 10"^ measurements of Wilson loops for each instanton configuration. Actually, we estimate 
the expectation value of Wilson loops from about 5 x 10^ measurements of Wilson loops on 
500 completely independent configurations. Figure |^ shows the Wilson loop for the case of 
z/ = 5. The Wilson loop seems to decay as {W[R,T]) oc exp[—aRT], which indicates the 
existence of a linear confining potential. As shown in Fig.^ the static quark potential V{R) 
is proportional to the inter-quark distance R up to the intermediate region R ~ 1.2 fm (< T). 
Here, we have checked that the static potential at i? < 1.2 fm does not depend on T (^ R) 
within errors. We speculate that instantons would be relevant degrees of freedom for the 
linear potential between a quark and anti-quark pair in the physically interesting region. 

The area low behavior of the Wilson loop means that the errors of this value grow as 
the area RT increases, since the Wilson loop decreases exponentially with RT. The static 
quark potential is not yet demonstrated at longer distance region as R > 1.2 fm, where it 
is necessary to calculate larger Wilson loops with T ^ 2.0 fm. To this end, we need a huge 
number of independent measurements and perform the smearing procedure to reduce the 
statistical error because of the small value of Wilson loops {W[R,T]) oc exp[— ct-RT]. 

Finally, we discuss the relation between the strength of color confinement and the in- 
stanton density. In order to extract the string tension, we measure the Creutz ratio 

(prj.. {WiR,T)){WiR-l,T-l)) 
^^ ' '~ {W{R,T-1)){W{R~1,T)) ^^ 

from the SU(2) Wilson loop. If the logarithm of the Wilson loop can be approximated as 
— ln{W{R, T)) = aRT + m{R + T) + const, the Creutz ratio gives the exponent of the string 
tension, x = exp(— cr^ai). The dimensionless lattice string tension aiat provides the physical 
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string tension aphys with the lattice-spacing a as Cphys = criat/ci^- 

We estimate the string tension around R ~ 0.8 fm for each instanton density, (N/V)* = 
0.50, 0.75 and 1.0 fm" . As shown in Fig.^ the string tension depends directly on the 
instanton density, and grows drastically as the instanton density increases. At a density of 
(N/V) = (1/fni)^ with z/ = 5, the string tension conies out to be about a ~ 0.4 GeV/fm, 
which corresponds to about a half of the physical string tension (Jexp — 0.89GeV/fm. 

In addition, we demonstrate the u = 3 case. As shown in Fig.^ and Fig.|lT], the Wilson 
loop seems to obey the area low and there appear a linear potential at R < 1.2 fm, which 



are qualitatively similar to the z/ = 5 case. Figure |T2| shows that the string tension of the 
u = 3 case is larger than that of the u = 5 case at the same density, which is caused by the 
large overlapping of instantons and anti-instantons in the u = 3 case. Rigorously speaking, 
for z/ = 3 it may not be suitable to adopt the simple sum ansatz, which is based on the 
statistical independence of the instanton ensemble. Therefore, one needs to construct a new 
formulation, which is workable also for the highly overlapping instanton system as the v = 3 
case. 

Finally, we would like to discuss the difference of the present results with the previous 
works. In our previous calculation [^], we fixed the peak of the instanton size distribution 
as Ppeak = 0.4 fm, which corresponds to the average instanton size of p ~ 0.45 fm. We 
extracted then the static potential using smaller Wilson loops. Since the string tension 
strongly depends also on the average instanton size p, we obtained a larger value for the 
string tension [^ than that of the present calculation. Recently, Brower et al. report that 



the slope of the heavy quark potential is about 0.1 GeV/fm at the average instanton size 
p = 1/3 fm and the density {N/V) = (l/fni)^ |^B[. The main reason of the discrepancy 
is the different choice of the average instanton size p. In addition, they use larger Wilson 
loops with T above 3 fm. The large Wilson loop justifies the extraction of the string tension. 
In this sense, they should obtain a better value for the string tension at long distance []68 



However, in their calculation, the perpendicular 2-dimensional directions {x — y) to the plane 
{z — t) of Wilson loops are limited to be much narrower than the size of the Wilson loop, 
although the other directions are very wide for the calculation of large Wilson loops. Our 
preliminary study of the contributions of instantons away from the plane of the Wilson loop 
in the x — y plane indicates appreciable contributions to such a large Wilson loop. Hence, 
it would be necessary to consider a large volume even in the x — y direction for the large 
Wilson loop, which is outside of the scope of this work. 

VII. SUMMARY AND CONCLUSION 

We have studied color confinement and nonperturbative quantities of the QCD vacuum 
using the multi-instanton configuration. We have made the present study by being motivated 
by the presence of a strong correlation between instantons and monopoles after the abelian 
gauge fixing in pure SU(2) gauge theory. In our calculation, the multi-instanton system is 
constructed by assuming the suppression of the large size instanton as /(p) ~ p^^ due to 
the infrared repulsive interaction between instantons. 

First, we have investigated the monopole-loop distribution in the multi-instanton system 
by using the maximally abelian gauge. Here, we have executed the "block-spin" transfor- 
mation in order to extract the infrared properties of monopole loops. The dilute-instanton 
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system produces small monopole loops localized around each isolated instanton. If instan- 
tons are well separated each other, monopole loops disappear in the continuum limit a — >■ 0. 
However, as the instanton density becomes higher, several small monopole loops combine 
into one longer loop. At a high instanton density, there appears a highly complicated 
monopole loop covering the entire physical volume. The appearance of long monopole loops 
of this system resembles that of the low-temperature lattice QCD where color confinement 



is realized through monopole condensation [48|. Our results indicate that instantons play 



an essential role on the promotion of a global network of monopoles. 

We have found that the high instanton density provides a highly complicated and long 
monopole loop. This seems to indicate that instantons are responsible also for confinement. 
Therefore, we have calculated the Wilson loop in the multi-instanton system without the 
abelian gauge fixing. Here, we take the instanton density of {N/V) = (1/fm)^ and the 
average instanton size of p = 0.4 fm. We have found that the instanton ensemble gives 
an area law behavior of the Wilson loop and the static quark potential is approximately 
proportional to the inter-quark distance, V{R) ~ aR, up to i? ^ 1.2 fm. In this situation, 
the string tension has been evaluated as cr ~ 0.4 GeV/fm, which corresponds to about a 
half of the physical string tension aexp — 0.89GeV/fm. Furthermore, we have discussed 
the dependence of the string tension on the instanton density. The string tension tends 
to decrease monotonously as the instanton density becomes smaller. Such a tendency is 
consistent with the disappearance of long monopole loop as the instanton number decreases. 
At the high temperature QCD vacuum where the string tension is zero, instantons and anti- 
instantons almost disappear and also the long monopole loops do. From this close relation 
between the monopole clustering and the instanton density, we speculate that instantons 
would be relevant degrees of freedom for the linear potential between quark and anti-quark 
in the physically interesting region. 

In the multi-instanton system, the string tension depends directly on both the instanton 
density and the average instanton size . Therefore, one has to fix these parameters from the 



original QCD vacuum. To this end, for instance, the improved cooling [|T2[ and the inverse- 
blocking |13|] are very interesting scheme for the extraction of topological quantities like 
instantons. Finally, it is necessary to simulate the static quark potential at longer distance 
(-R ^ 1.2 fm) using larger Wilson loops and enormous configurations in order to clarify the 
confinement properties up to very long distance. However, the calculation of such a large 
Wilson loop would be facing the limit of the present computation power. 
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FIGURE CAPTIONS 

FIG. 1. The instanton size distribution f{p) as a function of the instanton size p. The instanton 
size distribution for i/ = 5 is denoted by sohd curve and that for i/ = 3 by dashed curve. In both 
cases, the average size p is kept to p = 0.4 fm. 

FIG. 2. A schematic demonstration of the periodic condition used for an instanton ensemble 
sHced into 2-dimensional plane. We introduce the lattice in the center of the instanton ensemble 
to calculate the link variables by considering all the instantons outside of this central box obtained 
with a periodic condition for the sizes, positions color orientations of instantons. 

FIG. 3. The local correlation between instantons and monopoles for a dilute instanton system 
(a), and for a dense instanton system (b). The thick lines denote monopole currents which appear 
at some time t. The bottom surface corresponds to the action density of instantons s{x,y) at a 
certain time slice. Here, we show monopole currents which are transformed by one "block-spin" 
on the 32^ x 8^ lattice with a = 0.125 fm. 

FIG. 4. The histograms of monopole-loop length with various densities {N/VY'^ = 0.75, 1.0 
and 1.25 fm~ , after the various "block-spin" transformations in the v = 5 case, (a) denote the 
results of no "block-spin" transformation, (b), (c) and (d) denote the results after one, two and 
three "block-spin" transformations, respectively. 

FIG. 5. The histograms of monopole-loop length in unit of the lattice spacing a at high tem- 
perature (a), and at low temperature (b) in finite temperature SU{2) lattice QCD. 

FIG. 6. A random instanton configuration without the periodic condition for sizes, positions 
and color orientations of instantons sliced into 2-dimensional plane. We introduce the lattice in 
the center of the random instanton ensemble, a denotes the lattice unit and 36 corresponds to the 
number of lattice point in each time-space direction. 

FIG. 7. The ensemble average of the Wilson loop in the multi-instanton system with the infrared 
instanton size distribution of z^ = 5 as a function of the area R x T. The instanton density is 
(N/V)^^'^ = 1 fm~^ and the average size is p = 0.4 fm. 

FIG. 8. The static potential in the multi-instanton system as a function of the distance R 
extracted from the Wilson loop with T = 1.8 fm of FigJ^. 

FIG. 9. The string tension cr as a function of the instanton density (N/V)^'^ with the infrared 
instanton size distribution u = 5. 

FIG. 10. The ensemble average of the Wilson loop in the multi-instanton system with the 
infrared instanton size distribution z^ = 3 as a function of the area R x T. The instanton density 
and the average size are kept same as the case of z^ = 5. 

FIG. 11. The static potential in the multi-instanton system as a function of the distance R 
extracted from the Wilson loop of Fig.p!o|. 
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FIG. 12. The string tension o" as a function of the instanton density {N/VY''^ with the infrared 
instanton size distribution u = 2>. 
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